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1 Introduction 



A non-zero value of different QCD condensates leads to non-perturbative power correc- 
tions to propagators. The one being intensively studied during last years is the A^- 
condensate in Landau gauge [TJ |21 El IH El E] (extended to a gauge-invariant non-local 
operator, [7]), that is responsible for ~ 1/g^ corrections to the gluonic propagator com- 
pared to perturbation theory. In this paper we investigate the role of such corrections in 
the ghost propagator, and present a method that allows to test numerically that power 
corrections of ~ type really exist using only ghost and gluon lattice propagators, and 
ordinary perturbation theory. 

The study of the asymptotic behaviour of the ghost propagator in Landau gauge in 
the SU (3) quenched lattice gauge theory with Wilson action was the object of a previous 
work |H^. The lattice definition and the algorithm for the inversion of the Faddeev-Popov 
operator, as well as the procedure of eliminating specific lattice artifacts, are exposed 
there. A perturbative analysis, up to four- loop order (0,^01); has been accomplished 
over the whole available momentum window [2GeV ^ 6GeV]. However, a lesson we 
retained after a careful study of the gluon propagator performed in the past [HI [121 121 is 
that non-perturbative low-order power corrections and high-order perturbative logarithms 
give comparable contributions over momentum windows of such a width. Both appear to 
be hardly distinguishable, and thus - because of the narrowness of the fit window - the 
power-correction contribution could lead to some enhancement of the Aqcd parameter. 
Conversely, higher perturbative orders could borrow something to the non-perturbative 
condensate fitted from the power correction term. So, the quality of the fits (the value 
of x^/d.o.f) of lattice data is not a sufficient criterion when interpreting the results. A 
solution to the problem is to use several lattice data samples in order to increase the 
number of points in the fit window. This presumably brings another bias: the rescaling 
of the lattice data from different simulations (with different values of the ultraviolet (UV) 
cut-off i.e. the lattice spacing a). Nevertheless, we have to assume anyhow that the 
dependence on UV cut-off approximatively factorises ^ in order to fit lattice data to any 
continuum formula. Such an assumption will be furthermore under control provided that, 
as it happens in practice, our lattice data from different simulations match each other 
after rescaling. 

In the present paper we will follow the approach presented in refs. |2l [HI and do a fully 
consistent analysis of ghost and gluon propagators in the pure Yang-Mills theory based 
on the OPE description of the non-perturbative power corrections in Landau gauge. As 
far as our lattice correlation functions are computed in Landau gauge, the leading non- 
perturbative contribution is expected to be attached to the v.e.v. of the local oper- 
ator. This condensate generates a 1/g^-correcting term, still sizeable for our considered 
momenta, and that, as will be seen, gives identical power corrections to both gluon and 
ghost propagators. This result allows to separate the dominant power-correction term from 
the perturbative contribution, and suggests a new strategy for analysing the asymptotic 
behaviour of ghost and gluon propagators, even in the case of a small fit window. 

In the present letter we use this strategy to extract the Aqco-parameter from ghost 
and gluon propagators. 

'^This is the case of any renormalisation scheme where one drops any regular term depending on the 
cut-off away from renormahsation constants |13| . 
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2 The analytical inputs 

The present section is devoted to briefly overview the analytical (perturbative and non- 
perturbative) tools we have implemented to analyse our gluon and ghost lattice propaga- 
tors. 



2.1 Pure perturbation theory 

In the so-called Momentum subtraction (MOM) schemes, the renormalisation conditions 
are defined by setting some of the two- and three-point functions to their tree-level values 
at the renormalisation point. Then, in Landau gauge. 



dln(Z3,MOM(p^ = /i^ A) 
A^oo d In /i' 



lim -JTr—2 = 73,mom(5'mom) (1) 



where A is some regularisation parameter (a ^ if we specialise to lattice regularisation) 
and ^ 

^3.MOm(/ = ^^ A) = ■ P' ■ - ^) (-P)^^(P)) • (2) 

A similar expression can be written for the ghost propagator renormalisation factor Z3. 
Both anomalous dimensions for ghost and gluon propagators have been recently com- 
puted jH] in the MS scheme. At four- loop order we have 

dlniZs^MOAi) 13 , 3727 . /2127823 9747 



MS 



/3011547563 18987543 1431945 \, 
^ V 6912 256 64 ' 

ciln(Z3,MOM) _ 9 ^ ,813^2 , /^157303 5697 
4 

( 



dln'fl^^''^ - 4 Vs + 16 ^MS + 64 32 > '''^^ 



/219384137 9207729^ 221535 

+ C3 C5 ^ 

' 1536 512 32 ' 



(3) 



MS 



where h = g'^ /{AtxY. However, the definition of a MOM scheme still needs the definition 
of the MOM coupling constant. Once chosen a three-particle vertex, the polarisations 
and momenta of the particles at the subtraction point, there is a standard procedure to 
extract the vertex and to define the corresponding MOM coupling constant. This may 
be performed in several ways. In fact, infinitely many MOM schemes can be defined. In 
ref. the three- loop perturbative substraction of all the three- vertices appearing in the 
QCD Lagrangian for kinematical configurations with one vanishing momentum have been 
performed. In particular, the three schemes defined by the subtraction of the transversal 
part of the three-gluon vertex (MOMg) ^ and that of the ghost-gluon vertex with vanishing 
gluon momentum (MOMc) and vanishing incoming ghost momentum (MOMcO) will be 

"'in Euclidean space. 

■''It corresponds to MOMgg in |14| . 
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used in the following. In Landau gauge and in the pure Yang-Mills case {uf = 0) one has 

70 2 /51627 153 \ 3 

^MOM; -^MS + y /^MS + ( y ^3 I /iMS + 



/ 304676635 299961 81825 , ^ 
^ V~6912 6^^' " ~Qr^' ' ^ 



_ 223 2 /918819 351 > 3 

^MOM, ~ MS + ^"-MS + I ^296 ^ / 



/29551181 137199 74295 \ 
+ ^^64 32~^^ ~ ~6r 

169,2 /76063 153 \ 3 



''MS 



V'oSCo ~^MS + ^MS + ;l44 4 ^3 J ^MS+ 

/42074947 35385 66765 
+ ^^728 ~ ~6^^V ' 



(4) 



Thus, inverting Eq. (jH) and substituting in Eq. Q, we obtain the gluon and ghost propa- 
gator anomalous dimensions in the three above-mentioned renormalisation schemes. The 
knowledge of the /3-function 



7 n 

^^^^ = dhT^ h = -J2p, h^^^ + O (h-^^) , (5) 

makes possible the perturbative integration of the three equations obtained from Eq. (jH}. 
The integration and perturbative inversion of Eq. ^ at four-loop order gives an expression 
for the running coupling: 

1 A Alog(t) /9? 1 / / , , 1^ ^ 
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where t = In -7^ — . We omit the index specifying the renormalisation scheme both for h 

^^QCD 

and Aqcd- 

The last equation allows us to write the ghost and gluon propagators as functions of 
the momentum. The numerical coefficients for the /3-function in Eq. ^ are |15j : 

Po =}l^ Pi = 102, 

^MOM^ = 3040.48, /^a™" = 2412.16, pf^^ = 2952.73, 
^MoKo ^ 100541, pf^' = 84353.8, pf^^ = 101484. 

2.2 OPE power corrections for ghost and gluon propagators 

The dominant OPE power correction for the gluon propagator has been calculated in 
(120]), and it has the form 

Z,ie) - (1 + ^§^) . (7) 



In this section we present the calculation of the analogous correction to the ghost propa- 
gator. The leading power contribution to the ghost propagator 



F'^^q') = j (fxe"'-^ T [c\x)c\Q)'j ), (8) 

as in refs. [21 El for gluon two- and three-point Green functions, can be computed using 
the operator product expansion [1^] : 



T (nx)m) = J2 i^tt (x) 0,(0) 



(9) 



where Ot is a local operator, regular when x 0, and where the Wilson coefficient q 
contains the short-distance singularity. In fact, up to operators of dimension two, nothing 
but 1 and : : contribute to Eq. (jH)) in Landau gauge ^. Then, applying Q to (jH)), 

we obtain: 



pah I J2\ 



\ abcFT 



+ [c,]-' (g')(:^:(0)4(0):) 



pert V 



+ W 



ah 



+ 



(10) 



c 



where 



w 



ah 



/rfW"-^ (At'(O) T 



Km 



connected 



G(2)::,G(2) 



tt' 



2 X 



:iii 



and the SVZ sum rule ^\ is invoked to compute the Wilson coefficients. Thus, one should 
compute the ^^sunsef diagram in the last line of Eq. (fTT|l . that couples the ghost prop- 
agator to the gluon A^— condensate, to obtain the leading non-perturbative contribution 
(the first Wilson coefficient trivially gives the perturbative propagator). Finally, 



F''\q' 



Fpertil 



3 9UA')b. 



q- 



) + 0{g\q-') 



(12) 



where the A^-condensate is renormalised, according to the MOM scheme definition, by 
imposing the tree-level value to the Wilson coefficient at the renormalisation point, 
As far as we do not include the effects of the anomalous dimension of the operator (see 
ref. 0), we can factorise the perturbative ghost propagator. Then, doing the transverse 
projection, one obtains the following expression for the ghost dressing function: 



Z3{q' 



^3, pert 



3 9UA')r 
g2 4(iV2-r 



(13) 



^Those operators with an odd number of fields {d^A and d^c) cannot satisfy colour and Lorentz in- 
variance and do not contribute with a non-null non-perturbative expectation value, neither cc contributes 
because of the particular tensorial structure of the ghost-gluon vertex 
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We see that the multiphcative correction to the perturbative Z^ p^it is identical to that 
obtained in ref. [2j for the gluon propagator (Eq. ((Zj)). 

We do not know whether there is a deep reason for the equahty of the Wilson coeffi- 
cients at one loop for the gluon and ghost propagators. Is it a consequence of the absence 
of (gauge-dependant) {A^) contributions in gauge-invariant quantities? In principle, this 
could be proven either by a direct calculation of some gauge-invariant quantity or by 
analysing a Slavnov- Taylor identity ^H] that relates the ghost and gluon propagators with 
the three-gluon and the ghost-gluon vertices. In both cases one has to evaluate the (A^) 
corrections to these vertices, and this is a delicate question (because of soft external legs, 
[T^). The understanding of the mechanism of compensation of diverse gauge-dependent 
OPE contributions deserves a separate study, and we do not address this question in the 
present paper. 



3 Data Analysis 
3.1 Lattice setup 

The lattice data that we exploit in this letter were previously presented in ref. jH]. We 
refer to this work for all the details on the lattice simulation (algorithms, action, Faddeev- 
Popov operator inversion) and on the treatment of the lattice artifacts (extrapolation to 
the continuum limit, etc). 

The parameters of the whole set of simulations used are described in Tab. ^ 





Volume 


(GeV) 


Number of conf. 


6.0 


16^ 


1.96 


1000 


6.0 


24^ 


1.96 


500 


6.2 


24^ 


2.75 


500 


6.4 


32^ 


3.66 


250 



Table 1: Run parameters of the exploited data ([Hj)- 



Our strategy for the analysis will be, after rescaling and combining the data from each 
particular simulation, to try global fits over a momentum window as large as possible. As 
will be seen, after such a multiplicative rescaling, all the data match each other from ~ 2 
GeV to ~ 6 GeV (cf. Fig. Q). For the sake of completeness, we have furthermore per- 
formed an independent analysis (at fixed lattice spacing) for all simulations from Tab. ^ 
The results of this analysis are given in Appendix El 



3.2 Extracting Aqcd from lattice data 

Given that at the leading order the non-perturbative power corrections factorise as in 
Eq. (|T!^ and are identical in the case of the ghost and gluon j2| propagators, our strategy 
to extract Aqcd is to fit the ratio 

Zs{q^An),{A')) Z3,pert(g^A^)' ^ ' 
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to the ratio of three-loop perturbative formulae in scheme R obtained in section 12.11 and 
then convert Ar to Aj^ (|H1)- It is interesting to notice that non-perturbative correc- 
tions cancel out in this ratio even in the case nj ^ 0. The Aqco-parameter extracted 
from this ratio is free from non-perturbative power corrections up to operators of dimen- 
sion four, while the dressing functions themselves are corrected by the dimension two 
A^— condensate. In Tab. |21 the best-fit parameters for the three schemes are presented 
and we plot in Fig.Qthe lattice data and the MOM^ best-fit curve for the ratio in Eq. (fT^ . 



MOMg renormalizadon scheme 



MOMg renormalization scheme 




S 280 



k . (GeV) 



(b) 



Figure 1: (a) Plot of the |^ for the best fit parameter = 269(5) MeV. (b) The determi- 
nation of the optimal window fit (from 3 GeV to /cmaxO ^ """/S) results from the search for some 
"plateau" of Aj-jg when one changes the low bound of the fit window. 



scheme 


A 2 loops 


xVd.o.f 


A 3 loops 


xVd.o.f 


A 4 loops 


xVd-o.f 


MOMg 


324(6) 


0.33 


269(5) 


0.34 


282(6) 


0.34 


MOMc 


351(6) 


0.33 


273(5) 


0.34 


291(6) 


0.33 


MOMco 


385(7) 


0.33 


281(5) 


0.34 


298(6) 


0.33 



Table 2: The best-fitted values of Aj^ for the three considered renormalisation schemes. As 
discussed in the text, MOM^ seems to be the one showing the best asymptotic behaviour. 

In Fig. 121(a) we show the evolution of the fitted parameter A^^ when changing the 
order of the perturbation theory used in the fitting formula. One can conclude from 

Fig. 121(a) and App. ^that MOM^ scheme at three loops gives the most stable results for 
Aj^. It can also be seen from the ratio of four to three loops contributions (see Fig. 121(b)) 
for the perturbative expansion of logZ^, 

ln(Z3) = roln{hii) + Y,nh)i, (15) 

i=l 

where the coefficients are to be computed from those in Eqs. (jHliZI) and R stands for 
any renormalisation scheme {R = MOM^ in Fig. 121(b)). The same is done for logZ^. 
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According to our analysis, three loops seems to be the optimal order for asymptoticity. 
Indeed, the values of Aj^g for the three considered renormalisation schemes practically 
match each other at three loops. Finally, 

Ams = 269(5)^^2 (16) 

could be presented as the result for the fits of the ratio of dressing functions to perturbative 
formulae, where we take into account the bias due to the choice of the fitting window (see 
Fig.CKb) and App. Ej). 

However, there are indications (App. ^ and |8j) that our present systematic uncer- 
tainty may be underestimated, and we prefer simply to quote Aj^ !^ 270 MeV . This value 
is considerably smaller than the value of ~ 320 MeV obtained by independent fits of dress- 
ing functions (jH]), and with fitting windows independently determined for each lattice 
sample (see Fig. 121(a))- This argues in favour of presence of low-order non-perturbative 
corrections to the ghost and gluon propagators. 



number of loops 



(a) 



• ratio MOMg 
■ ratio MOMc 
' ratio MOMcO 

♦ Gluon MOMg 

* Ghost MOMg 



(r,h)/rj 




3 3.5 4 4.5 S 5.5 6 3 3.5 4 45 5 5.5 6 



k (GeV) 



k (GeV) 



(b) 



Figure 2: (a) Evolution of the parameter Aj^, extracted from fits of the ratio Eq. (|14j) and 
propagators alone (rhombus and star markers, extracted from Tab. 7 of [H]) to perturbative 
formulae, as function of the order of the perturbation theory. Only statistical error is quoted, 
(b) Ratio of four-loop to three-loop contributions (and of three-loop to two-loops for the sake of 
comparison) for the perturbative expansion of logZ^ and logZ^ (in MOM^) in Eq. ()15|). plotted 
versus the momenta inside our fitting window. 



3.3 Estimating the value of the (^4^) gluon condensate 

Knowing Aj^ we can fit ghost and gluon dressing functions using Eqs. ()7|13j) . The 
free parameter in this case is g\{A^). According to the theoretical argument given in 
12.21 the results obtained from these fits have to be compatible. We have performed this 
analysis for the rough value Kyi^ ~ 270 MeV, (see Tab. |HI). Indeed, we find that the 
resulting values agree. It is worth to emphasise the meaning of this result: a fully self- 
consistent description of gluon and ghost propagators computed from the same sample of 
lattice configuration (same A-^ and same {A^)) is obtained. 
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^3 


gliA') (GeV^) 


2.7(4) 


2.7(2) 



Table 3: The best-fitted values of gj^{A'^) for MOM^ obtained from fitting lattice data to a 
three-loop perturbative formula + non-perturbative power correction with = 270 MeV. We 
only quote statistical errors. 

The values of the gluon condensate presented in Tab. Elare smaller than those obtained 
from the previous analysis of the gluon propagator j2]. The reason for this is the larger 
value of we have found. Had we taken Aj^ ~ 240 MeV, we would obtain similar 
results to those previously presented. Of course, this discrepancy have to be included 
in the present systematical uncertainty of our analysis of the ghost propagator lattice 
data. However, the purpose of this paper is not to present a precise determination of 
the dimension-two gluon condensate, but only to show that ghost and gluon propagators 
analysis strongly indicates its existence. The precision could be improved by increasing 
the Monte-Carlo statistics and by performing new simulations at larger (3. 

Another source of discrepancy are renormalon-type contributions that can also be of 
the order of ~ In fact, our OPE study does not include the analysis of such correc- 

tions. However, the numerical equality (cf. TabOj) of ~ power corrections at fixed 
Ayf^ suggests that the ratio (fT^ is free of such corrections, in agreement with the com- 
mon belief that the renormalon ambiguities are compensated by condensate contributions. 
The estimate for Aj^ obtained from this ratio is thus not affected by the renormalon-type 
contributions. But the dependence of the value of these corrections on Aj^ speaks in 
favour of the presence of renormalon-type contributions in Z3 and Z3 separately. 

4 Conclusions 

We have analysed non-perturbative low-order power corrections to the ghost propagator 
in Landau gauge pure Yang-Mills theory using OPE. We found that these corrections are 
the same as those for the gluon propagator at leading order. This means that their ratio 
does not contain low-order power corrections (~ 1/?^), and can be described (up to terms 
of order ~ ^/q^) by the perturbation theory. Fitting the ratio of propagators calculated on 
the lattice we have extracted the Aj^ parameter using three- and four-loop perturbation 
theory. The value Aj^ ^ 270 MeV extracted from the ratio is quite small compared to 
the one obtained in fits of gluon and ghost propagator (A£^* ~ 320 MeV, |H]) separately. 
Indeed, Aj^ ^ 270 MeV extracted from the ratio of ghost and gluon dressing functions 
is closer to the value calculated in the past with power-corrections taken into account 
(A^ii^^ ^ 250 MeV, |S1I11) than to the purely perturbative result . This study within 
perturbation theory confirms the validity of our OPE analysis, and argues in favour of a 
non-zero value of non-perturbative A^-condensate. We are not able at the moment to give 
a precise value of the A^-condensate using this strategy. More lattice data and detailed 
analysis of diverse systematic uncertainties are needed for this. But the method exposed 
in this letter can in principle be used for this purpose, both in quenched and unquenched 
cases. 
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A Appendix 



In this Appendix we present results for Aqcd extracted by fitting the ratio using 
two, three and four-loop perturbation theory. But we do not mix data samples obtained 
in different lattice simulations. This allows to control the effects of several lattice artifacts 
and of the uncertainty on the lattice spacing calculation on the resulting value of Aqcd- 

Fits have been performed in MOMc, MOMco, MOM^ renormalisation schemes (cf. Tab. 
mHEEIHEl) • In each case we chose the best fit from several fitting windows, having the 
smallest x^/d-o.f.; the statistical error corresponds to that fit. The systematic error is 
calculated from different fit windows. 

One can see from these fits that the values for Aj^ are small at three and four loops 
when fitting at energies > 3GeV. All results are rather stable in this domain, and thus the 
fitting of combined data from the simulations with different lattice spacings, presented in 
the main part of the present letter, is safe and well defined. 
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Fits at two loops 



V 




Left, GeV 


Right, GeV 


aAy 

|a,MOMg 


conversion to A^"'^^^^^, MeV 


x2/d.o.f. 




6.0 


2.54 


4.32 


529(17)+^ 


359(12)1^ 


0.21 




6.0 


3.14 


4.12 


513(15)+^^ 


348(10)1^1 


0.10 


24^ 


6.2 


3.02 


4.95 


377(24)_n 


358(22)_io 


0.14 


32^ 


6.4 


3.66 


5.85 


257(21)+^ 


325(26)1^ 


0.17 



Table 4: Perturbative fit of r, ,1 at 2 loops in MOMo scheme and further conversion to 



V 




Left, GeV 


Right, GeV 


||,MOMc 


conversion to A(2)ms^ jyjgy 


x2/d.o.f. 


16^ 


6.0 


2.15 


4.12 


445(6)_6 


375(5)_5 


0.14 


24"^ 


6.0 


3.14 


4.12 


398(53)li^ 


335(45)li' 


0.10 


24^ 


6.2 


3.02 


4.95 


313(19)_22 


369(22)„26 


0.13 


32^ 


6.4 


3.66 


5.85 


215(17)1^ 


337(26)1;] 


0.17 



Table 5: Perturbative fit of Zz(p'2) ^ loops in MOMc scheme and further conversion to 
MS 



V 




Left, GeV 


Right, GeV 


aAg^ ^ 

||,MOMco 


conversion to A^^^'^^, MeV 


x2/d.o.f. 


16^ 


6.0 


1.97 


4.11 


400(6)„5 


413(6)_5 


0.15 


24"^ 


6.0 


3.13 


4.12 


354(49)+26 


367(41)+27 


0.11 


24^ 


6.2 


3.02 


4.95 


280(17)1^2 


367(24)1^7 


0.11 


32"^ 


6.4 


3.66 


5.85 


190(16)1^ 


366(30)1^ 


0.16 



Table 6: Perturbative fit of at 2 loops in MOM^o scheme and further conversion to 

MS 



11 



Fits at three loops 



V 


/5 


Left, GeV 


Right, GeV 


A (3) 
|i,MOMc 


conversion to A^"*^^^^*, MeV 


9/1 r 

X /d.o.f. 




6.0 


2.54 


4.31 


354(12)+^ 


297(10)+! 


0.23 




6.0 


3.13 


4.12 


312(48)+^° 


261(40)+25 


0.10 


244 


6.2 


3.14 


4.95 


247(20)„22 


289(23)^26 


0.14 


26 32^ 


6.4 


3.66 


5.86 


163(15)+^ 


254(24)+^ 


0.16 



Table 7: Perturbative fit of z-i(p^) ^ loops in MOMg scheme and further conversion to 
MS 



V 




Left, GeV 


Right, GeV 


||,MOMg 


conversion to A^^)'^^, MeV 


xVd.o.f. 


16^ 


6.0 


2.06 


3.44 


453(18)_5 


307(12)_4 


0.16 


24"^ 


6.0 


2.06 


3.53 


451(16)^^0 


306(11)+^ 


0.1 


24^ 


6.2 


2.80 


4.81 


295(20)+^^ 


282(19)+^^ 


0.13 


32^ 


6.4 


3.79 


5.59 


216(25)_i3 


270(31)„i6 


0.13 



Table 8: Perturbative fit of ^ loops in MOM^ scheme and further conversion to 

MS 



V 




Left, GeV 


Right, GeV 


aA?^ ^ 

||,MOMco 


conversion to A^^^'^^, MeV 


xVd.o.f. 


16^ 


6.0 


2.16 


3.53 


312(ll)_i7 


323(ll)_i8 


0.11 


24"^ 


6.0 


3.21 


4.11 


252(47)+36 


260(48)+37 


0.10 


24^ 


6.2 


3.13 


4.95 


205(16)_i8 


297(23)^26 


0.14 


32"^ 


6.4 


3.66 


5.86 


136(12)+^ 


262(23)+^ 


0.17 



Table 9: Perturbative fit of 7^7^ at 3 loops in MOM^o scheme and further conversion to 
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Fits at four loops 



V 




Left, GeV 


Right, GeV 


A (4) 

aAy 

|a,MOMg 

i_ 


conversion to A^^J^b^ j^^y 


x2/d.o.f. 




6.0 














6.0 


3.14 


4.12 


365(13)+!^ 


248(8)+!' 


0.10 


24^ 


6.2 


3.02 


4.95 


288(17)+^ 


274(16)+^ 


0.13 


32^ 


6.4 


3.66 


5.85 


199(15)+^ 


252(19)+^ 


0.17 



Table 10: Perturbative fit of Zzip"!) ^ loops in MOM^ scheme and further conversion to 
MS 



V 




Left, GeV 


Right, GeV 


aKf ^ 

||,MOMc 


conversion to A^^)'^^^ jy^^y 


x2/d.o.f. 


16^ 


6.0 












24"^ 


6.0 












24^ 


6.2 












32^ 


6.4 


3.66 


5.85 


175(15)+^ 


274(23)+^ 


0.17 



Table 11: Perturbative fit of Zz{p'2) ^ loops in MOMc scheme and further conversion to 
MS 



V 




Left, GeV 


Right, GeV 


aKf ^ 

yi,MOMco 


conversion to A^^^"^^, MeV 


x2/d.o.f. 


16^ 


6.0 












24"^ 


6.0 


2.95 


4.12 


299(29)_26 


309(30)_27 


0.11 


24^ 


6.2 


3.02 


4.95 


225(14)+^ 


326(20)1^ 


0.13 


32"^ 


6.4 


3.66 


5.85 


152(13)+^ 


293(25)+^ 


0.16 



Table 12: Perturbative fit of at 4 loops in MOM^o scheme and further conversion 

to MS 
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